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Abstract 

Gas-surface  vibrational  energy  accommodation  in  a  reactive 
system  has  been  studied  under  single  collision  conditions.  The 
reaction  system  is  the  isomerization  of  cyclobutene  to  1,3-butadiene 
Both  seasoned  pyrex  and  silica  surfaces  were  used  over  the  tempera¬ 
ture  ranges  400-775  K  and  500-975  K,  respecti vely .  Strong  collider 
behavior  was  observed  below  425  K.  The  vibrational  energy  accommo 
dation  coefficient  was  found  to  decrease  with  rise  of  temperature 
from  o,  1.0  to  0.2,  while  the  relative  collisional  efficiency 
declined  from  %  1.0  to  0.008.  A  stochastic  model  was  used  to  fit 
the  data. 
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Introduction 

Gas-surface  vibrational  energy  accommodation  has  been  a  subject  of 

1  ? 

considerable  interest  both  experimentally  and  theoretically.  *  However,  to 

this  time  there  have  been  no  single-collision  measurements  of  gas-surface 

3 

energy  transfer  in  a  high  temperature  reactive  system.  In  the  present  communi- 

4-8 

cation,  we  describe  the  application  of  the  Variable  Encounter  Method  (VEM)  to 
the  study  of  single  collision  measurements  with  seasoned  heated  glass  surfaces. 
Although  such  a  surface  is  not  well  defined,  it  is  the  conventional  experimental 
surface  of  thermal  unimolecular  kinetics. 

The  VEM  technique  has  previously  been  applied  to  the  study  of  energy  tran¬ 
sients:  molecules  are  thermally  equilibrated  at  some  low  temperature  before 
experiencing  an  encounter,  i.e.  a  sequence  of  m  collisions  (where  m  is  variable 
and  is  chosen  to  be  1  in  the  present  study)  with  a  hot  reactor  surface;  after 
the  encounter,  unreacted  molecules  are  reequilibrated  at  the  original  low  tempera¬ 
ture.  The  amount  of  gas  phase  unimolecular  reaction  that  occurs  following 
surface  activation  gives  the  fraction  of  those  molecules  in  the  original 
Boltzmann  distribution  Nc  that  have  been  transported  to  energy  levels  above 
the  critical  reaction  threshold,  Eq.  As  m  increases,  Nc  approaches  N^,  the 
steady  state  population  distribution  characteristic  of  reactive  depletion  at 
the  reactor  temperature. 

In  the  present  work,  the  hot  surface  is  convex  in  shape,  ensuring  that 
only  a  single  collision  occurs  on  encounter  of  the  molecule  with  the  surface. 

The  single  collision  case  has  the  advantage  over  previous  mul ticoll ision  VEM 

oo 

work  (where  m  =  \  nf(n)),  in  that  no  averaging  over  some  normalized  distri- 

n=l 

bution  f(n)  of  numbers  of  collisions  made  by  the  molecule  per  encounter  with 


the  reactor  is  necessary  in  the  deconvolution  of  the  data.  This  provides 
greater  experimental  stringency  and  constraints  on  suitable  models  for 
the  collisional  transition  probability  matrix,  P  .  The  present  study  also 


2 


permits  a  more  direct  evaluation  of  the  gas-surface  vibrational  energy  accommo¬ 
dation  coefficient,  a  . 

The  ring-opening  unimolecular  isomerization  of  cyclobutene  to  1,3  buta- 
9  «i 

diene  has  a  low  Eq  (32.4  kcal  mole  )  and  is  substantially  in  the  low  pressure 

-4 

region  at  typical  experimental  pressures,  1  -  5x10  torr.  Toe  reaction  pro- 

7,8 

ceeds  with  no  indication  of  surface  catalysis.  On  the  basis  of  the  increasing 
strength  of  collision  with  decreasing  temperature  that  was  observed  in  earlier 
VEM  studies,  these  characteristics  make  this  species  suitable  for  single 
collision  study,  and  a  significant  amount  of  reaction  should  occur  at  lower 
temperatures . 
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Experimental 

The  reactant  cyclobutene  was  99.9%  pure.  The  reactor  is  shown  schemati¬ 
cally  in  Fig.  1.  The  reaction  vessel  consisted  of  a  3  l  pyrex  flask  with  an 
internally  heated,  central  reactor  finger,  A.  Two  different  reaction  vessels 
were  used,  fitted  with  pyrex  and  silica  reaction  fingers,  respectively.  The 
vessel  wall  was  cooled  by  an  ice  bath.  Liquid  Bi-Pb-Sn  eutectic  was  used  as 
the  internal  bath  medium  (D)in  the  reactor  finger  and  the  top  and  bottom  regions 
were  heated  independently  (C)  to  minimize  axial  temperature  variations.  At 
the  highest  temperature  (970  K)  the  maximum  deviation  was  ±  5°C  and  was  signi¬ 
ficantly  less  at  lower  temperatures.  Temperature  measurements  were  made  with 
calibrated  chromel-alumel  thermocouples  (B).  Before  kinetic  measurements  were 
made,  the  reactor  was  "seasoned"  at  750  K  by  pre-exposure  for  several  hours 

_3 

to  cyclobutene  at  a  pressure  of  n,  5x10  torr. 

The  system  was  run  in  both  static  and  flow  modes.  In  the  former,  the 

reactor  was  pumped  to  <  10”6  torr  and  substrate  was  admitted  at  a  pressure 
-4  -3 

between  10  -  10  torr.  The  gas  was  pumped  out  and  analyzed  after  a  pre¬ 

determined  run  time.  For  flow-mode  runs,  the  impedance  of  a  variable  exit 

control  valve  was  calibrated  by  pressure-drop  measurements.  A  substrate  entry 

-2  1 

flow  control  provided  gas  delivery  rates  of  2.5x10  to  0.22  torr  cc  s  . 

Typical  residence  times  were  10  to  100  s,  and  pressure  rose  from  <  10  6  to 

_  3 

<  10  torr  during  a  run,  which  was  terminated  by  closing  the  variable  flow 
stopcock  after  about  30  to  60  s  .  Pressure  measurements  were  made  with  an 
MKS  Model  170  capacitance  manometer. 

Product  analysis  was  performed  by  gas-liquid  chromatography  on  a 
5'  x  3/16"  squalane  column  at  0°C  with  FID. 
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Results  and  Discussion 

The  probability  of  reaction  per  collision,  Pc  ,  was  calculated  from  the 

product  yield  in  both  static  and  flow  (Appendix)  cases,  and  is  plotted  as  a 

function  of  temperature  in  Fig.  2.  Uncertainties  in  Pc  varied  from  about  ±  7« 

at  the  highest  to  ±  20%  at  the  lowest  temperatures.  This  is  due  to  the 

greater  temperature  sensitivity  at  the  lower  temperatures.  Also  shown  is  a  strong 

collider  curve,  calculated  on  the  basis  of  the  vibration  frequency  assignments  of 
g 

Elliott  and  Frey.  The  experimental  curves  coalesce  at  lower  temperature  where 
strong  collider  behavior  is  approached.  At  higher  temperatures,  slightly  larger 
reaction  probabilities  were  obtained  with  the  pyrex  reactor  finger.  Although  the 
difference  between  the  experimental  curves  is  small,  it  exceeds  the  experimental 
uncertainty  at  the  higher  temperatures,  and  appears  to  represent  either  a  minor 
difference  in  the  rate  of  energy  transfer  on  the  pyrex-based  vy  silica-based 
surfaces  or,  possibly,  a  small  difference  in  the  preparation  and  seasoning  pro¬ 
cedures  for  the  two  fingers. 

The  mathematical  models  usually  used  to  interpret  energy  transfer  data 

10  4 

both  in  conventional  thermal  and  VEM  systems  give  the  probability.  P(AE),  of 
a  down-transition  of  magnitude  AE  by  the  substrate  molecule  as  an  exponential 
(eq.  1)  or  gaussian  (eq.  2)  function: 

P(AE)  =  c-j  exp  (-AE/<AE>)  (1) 

or 

P ( AE )  =  c2  exp  (-(AE-AEmp)2/2o2)  (2) 

where  <AE>,  AEmp  and  o  are  adjustable  parameters  and  the  c's  are  constants. 

These  functions,  along  with  the  detailed  balance  and  completeness  relations, 
specify  P.^  However,  these  forms  are  inadequate  to  describe  near-strong 
collider  behavior.  As  strong  collider  behavior  is  approached,  the  columns  of 
P  should  approach  normalized  Bol  tzmann  vectors  B  characteristic  of  the  reactor 
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tempera ture,Tr,  i.e.  p.^  =  B^,  for  all  i,  j;  this  condition  is  clearly  not 
met  by  eq.  1  or  2.  All  choices  of  parameters  with  these  models  give  values 
of  Pc  that  are  significantly  below  strong  collider  and  the  present  observed 
values.  The  Boltzmann-weighted  exponential  (BE)(eq.  3)  does  fit  the  data 
and  gives  the  correct  strong  collider  limiting  behavior;  here 

PlVj  =  c3  exp  (-AE  /<AE>)  gi  exp  (-E./kTr)  (3) 

g^  is  the  density  of  vibrational  states  at  level  E^ , 
and  c^  is  a  constant. 

The  calculated  average  amount  of  energy  <AE‘>  transferred  in  a  down 

E 

O 

step  from  hot  molecules  at  the  level  of  energy  Eq  is  shown  in  Table  II  for 
several  reactor  temperatures.  Still  another  quantity  of  interest  is  the 
average  energy  transfer,  AEaw,  defined  as  the  difference  in  average  energies 
of  substrate  molecules  before,  E. ,  and  after, Ef,  collision;  Ei  refers  to  the 
flask  temperature  (273K).  This  is  related  to  the  vibrational  accommodation 
coefficient  a  =  AE  /(E  -E-),  where  E  is  the  average  energy  of  substrate 
molecules  at  the  reactor  finger  temperature.  Values  of  a  and  of  AEflv  are  given 
in  Table  I.  Also  given  are  values  of  the  collisional  efficiency  B-j ,  defined  as 
$i  =  Pc/P(sc),  as  well  as  efficiencies  B$s  that  would  be  observed  in  an  equiva¬ 
lent  steady  state  thermal  system;  B$s  =  P(°°)/P(sc) ,  where  P(°°)  is  the  actual 
steady-state  reaction  probability  per  collision  and  P(sc)  is  the  strong-collider 
value  of  the  reaction  probability  per  collision.  Bss  is  the  same  value  that 
would  apply  in  a  low-pressure  gas-phase  thermal  unimolecular  system  where  it  is 

defined  as  8  =  k  wc/k  sc;  ojWC  =  wsc;  or  equivalently  as  B  =  uSC/u)WC  , 

ss  o  o  J  ss 

wc  sc  wc  sc 

k  =  k  ;  here  k  ,  k  are  the  observed  low  pressure  rate  constants 
oooo 

WC  s  c 

measured  at  constant  collision  frequencies  u>  and  u>  in  a  weak  collider  and 
strong  collider  system,  respectively.  It  is  evident  from  Table  I  that  since 
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the  Ss$  values  group  close  to  unity  at  all  temperatures,  a  less  stringent  test 

of  energy  transfer  behavior  is  available  in  steady-state  experiments  as  compared 

to  the  present  measurements  of  B-j  in  the  transient  region.  In  reactive  systems, 

B-|  is  also  a  more  sensitive  function  of  energy  transfer  efficiency  than  is  a  , 

since  the  former  tests  efficiency  relative  to  the  difference,  Eq-E^  ,  while  the 

latter  tests  efficiency  relative  to  the  difference,  E  -  E^ . 

The  decline  in  collisional  efficiency  that  occurs  with  rise  of  temperature 

7  8 

in  Fig.  2  and  Table  I  has  been  observed  previously  in  cyclobutene  VEM  studies  ’ 
where  strong  collider  behavior  was  inferred  below  %  450  K,  as  has  been  found 
here.  This  same  general  trend  has  been  observed  in  all  VEM  studies  and  in  con¬ 
ventional  homogeneous  thermal  studies  at  higher  temperatures.  It  is  also 
generally  observed  in  energy  accommodation  studies.^  This  trend  appears  to  be 
related  to  the  fact  that  as  the  temperature  is  lowered  the  time  that  the  mole¬ 
cule  spends  trapped  in  the  gas-surface  potential  well  increases.  The  Boltzmann- 
weighted  models  then  correspond  more  closely  to  the  case  of  statistical  accommo¬ 
dation  on  the  surface, with  relaxation  of  the  restrictions  on  conservation  of 

1 2 

angular  momentum  that  apply  in  gas-gas  systems. 

Draper  and  Rosenblatt  measured  vibrational  energy  accommodation  coefficients 

of  hydrocarbons  on  metal  surfaces  at  room  temperature  by  a  vibrating  diaphragm  tech- 
2 

mque.  Although  those  experiments  are  not  directly  comparable  to  ours,  we  note 

that  their  values  of  avib  are  fairly  close  to  unity,  i.e.,  0.7  to  0.8.  For 

larger  molecules,  i.e.,  as  the  depth  of  the  potential  well  increased,  they 

found  that  the  values  of  a  also  increased. 

vib 

One  final  comment:  the  present  data  on  energy  transfer  efficiency  extend 
upward  from  425  K;  the  low  temperature  region  from  %  200  to  %  700  K  is  one 
that  has  previously  been  explored  in  thermal,  chemical  activation  and  photo¬ 
excitation  reactive  systems  with  conflicting  conclusions  as  to  whether  effi¬ 
ciency  increases  or  decreases  with  rise  of  temperature.^  The  present  system 
reveals  marked  decrease  in  <AE'>,  and  a  in  this  region  with  rise  of  temperature. 
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Table  I.  Energy  transfer  parameters  in  the  cyclobutene  system 

T(K) 


400 

500 

600 

750 

875 

975 

<AE1>e  fcm"1) 

o 

*  - - 

10,550 

8940 

6925 

(7520) 

4200 

(4700) 

3125 

2540 

AE  (ciiT^) 
av  ' 

465 

720 

780 

(900) 

900 

(1000) 

1050 

1135 

a 

1.0 

0.57 

0.40 

(0.47) 

0.28 

(0.31) 

0.23 

0.20 

6SS 

1.0 

1  .0 

0.994 

(0.998) 

0.87 

(0.92) 

0.68 

0.51 

61 

1.0 

0.40 

0.14 

(0.22) 

0.029 

(0.048) 

0.0134 

0.008 

a)  Parenthetic  values  apply  to  the  pyrex  reactor.  Pyrex  and  silica  reactors 
give  indistinguishable  results  at  500  K. 

b)  Strong  collider  behavior 

i 

,1 
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Appendix.  Flow  system:  apparent  rate  constant. 

Mass  conservation  requires 

V  ^  =  Q‘  -  F'R  -  kVR  and  V  =  _F'P  +  kVR, 
dt  dt 

where  R(t)  and  P(t)  are  the  concentrations  of  reactants  and  products  at 
time  t,  respectively;  Q‘  is  the  flux  of  reactant  into  the  reaction  vessel 
of  volume  V;  k  is  the  apparent  rate  constant  (=  Pc  c  A/hV,  where  c  is  the 
mean  molecular  speed  and  A  is  the  heated  reactor  area);  and  F‘  is  the  Knudsen 
output  flow  impedance.  Since  the  reactor  in  initially  evacuated,  P(0)  =  0, 

R(0)  =  0.  The  amount  of  product  collected,  Py  ,  is  given  by 

PT  =  F '  P(t)dt 

0  t 

where  tffl  is  the  run  time.  An  analogous  expression  applies  for  the  amount 
of  reactant  collected  Ry.  One  obtains 

RT  =  F  Sn  -  4k  n-exp(-(F^)tm)] 

Ry+Py  F+k  t  _  1  [1_eXp(.F  t  )] 
m  F  L  r  m,J 

where  F  =  F'/V.  This  equation  may  be  solved  iteratively  for  k. 
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Figure  1. 


Figure  2. 


Figure  Captions 

Schematic  of  single  collision  VEM  reactor.  Shown  are  pyrex 
(or  silica)  reactor  finger.  A;  moveable  chromel-alumel  thermo¬ 
couple,  B;  heaters,  C,  in  two  sections,  top  and  bottom;  eutectic 
bath,  0;  and  variable  flow  stopcock,  E. 

Plots  of  experimental  values  of  Pc  vs  T  for  the  following 
reactors  and  conditions:  pyrex,  static,  fl  ;  silica,  static,  Q 
and  silica,  flow,  Q  .  Also  shown  is  a  calculated  strong 
collider  (SC)  curve.  The  magnitudes  of  random  errors  are 
approximately  given  by  the  size  of  the  symbols.  The  magnitude 
of  unspecified  systematic  errors  may  be  judged  from  the  scatter 
of  the  data;  such  errors  are  largest  at  low  temperatures  and  are 
believed  to  be  primarily  due  to  error  in  temperature  (see  text). 
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were  used  over  the  ,temperatnfe"ranges  400-775/ k  and  50D~^v^.rF;4reipei.i.i*cij. 
Strong  coJi^terTSehavior  was  observed  below©  425  K.  The  vibr3tt««4lenergy 
accommodation  coefficient  was  found  to  decrease  with  rise  of  temp  eratill'a.  from 
©fl.O  to  0.2,  while  the  relative  collisional  efficiency  /Qdecl  ined  fromcJl.Q 
to  0.008.  A  stochastic  model  was  used  to  fit  the  data.^r'  i 
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